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Abstract 

Recent research in generalizing quantum computation from 2-valued qudits to d- 
valued qudits has shown practical advantages for scaling up a quantum computer. 
A further generalization leads to quantum computing with hybrid qudits where 
two or more qudits have different finite dimensions. Advantages of hybrid and d- 
valued gates (circuits) and their physical realizations have been studied in detail by 
Muthukrishnan and Stroud (Physical Review A, 052309, 2000), Daboul et al. (J. 
Phys. A: Math. Gen. 36 2525-2536, 2003), and Bartlett et al (Physical Review A, 
Vol.65, 052316, 2002). In both cases, a quantum computation is performed when 
a unitary evolution operator, acting as a quantum logic gate, transforms the state 
of qudits in a quantum system. Unitary operators can be represented by square 
unitary matrices. If the system consists of a single qudit, then Tilma et al (J. Phys. 
A: Math. Gen. 35 (2002) 10467-10501) have shown that the unitary evolution matrix 
(gate) can be synthesized in terms of its Euler angle parametrization. However, if 
the quantum system consists of multiple qudits, then a gate may be synthesized 
by matrix decomposition techniques such as QR factorization and the Cosine-sine 
Decomposition (CSD). In this article, we present a CSD based synthesis method for 
n qudit hybrid quantum gates, and as a consequence, derive a CSD based synthesis 
method for n qudit gates where all the qudits have the same dimension. 
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1 Introduction 



A qudit replaces a classical dit as an information unit in ci-valucd quantum 
computing. A qudit is represented as a unit vector in the state space, which is a 
complex projective d dimensional Hilbert space, Hd- In the computational ba- 
sis, the basis vectors of Hd are written in Dirac notation as |0) , |1) , . . . |d — 1), 
where \i) = (0, 0, . . . , 1, . . . , 0)^ with a 1 in the {i + l)-st coordinate, for 
< i < {d — 1). An arbitrary vector \a) in Tid can be expressed as a lin- 
ear combination \a) = J2iZo l^), Xi & C and J2 = 1. The real number 
\xif is the probability that the state vector \a) will be in i-th basis state upon 
measurement. 

When the state spaces of n qudits of different d- valued dimensions are com- 
bined via their algebraic tensor product, the result is a n qudit hybrid state 
space H = Tidi ® ^^^2 ® ■ ■ ■ ® "^dn , where Tid^ is the state space of the rfj-valued 
qudit. The computational basis for Ti would consist of all possible tensor prod- 
ucts of the computational basis vectors of the component state spaces l-id^ ■ If 
di — d for each i, the resulting state space Tif " is that of n d- valued qudits. 

The evolution of state space changes the state of the qudits via the action of 
a unitary operator on the qudits. A unitary operator can be represented by a 
unitary evolution matrix. For the hybrid state space Ti, an evolution matrix 
will have size {did2 ■ ■ ■ dj^) x {did2 ■ ■ ■ dj^f), while the evolution matrix for 'Hf" 
will be of size d"" x d". In the context of quantum logic synthesis, an evolution 
matrix is a quantum logic circuit that needs to be realized by a universal set 
of quantum logic gates. It is well established that sets of one and two qudit 
quantum gates are universal [3,5,10, 15]. Hence, the synthesis of an evolution 
matrix requires that the matrix be decomposed to the level of unitary matrices 
acting on one or two qudits. 

Unitary matrix decomposition methods like the QR factorization and the Co- 
sine Sine decomposition from matrix perturbation theory have been used for 
2-valued and 3-valued quantum logic synthesis. In these domains, qudits are 
referred to as qubits and qutrits respectively. The Cosine Sine decomposition 
(CSD) of a unitary matrix, discussed in section 2, has been used by Mottonen 
et. al [9] and Shende et. al [12] to iteratively synthesize multi-qubit quantum 
circuits. The authors of this article recently extended the CSD to iterated 
synthesis of 3-valued quantum logic circuits acting on n qutrits [8]. Bullock 
et.al have recently presented a synthesis method for n qudit quantum logic 
gates using a variation of the QR matrix factorization [4] . This article presents 
a CSD based method for synthesis of n qudit hybrid and d-valued quantum 
logic gates. 
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2 The Cosine-Sine Decomposition (CSD) 



Let the unitary matrix WE C"*^"* be partitioned in 2 x 2 block form as 



W 



r m — r 



m — r \W' 



21 



22 



(2.0.1) 



with 2r < m. Then there exist r x r unitary matrices U and X, r x r real 



diagonal matrices C and S, and 
such that 



(m — r) X (m — r) unitary matrices V and F 



W 









V 




(2.0.2) 
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The matrices C and are the so-called cosine- sine matrices and are of the 
form C = diag(cos^^i, cos^^2! • • • • cos^^.), S = diag(sin^i, sin ^2, • • • , sin^^) such 
that sin^ 9i + cos^ 9i = 1 for some 9i, 1 <i <r [13]. Algorithms for computing 
the CSD and the angles 6i are given in [2, 14]. The CSD is essentially the 
well known singular value decomposition of a unitary matrix implemented 
at the block matrix level [11]. In sections 3 and 4, we give an overview of 
the CSD based synthesis methods of 2 and 3-valucd quantum logic circuits, 
respectively. From now on, we will not distinguish between gates, circuits and 
their corresponding unitary matrices. 



3 Synthesis of 2-valued (binEiry) Quantum Logic Circuits 



As shown in [8,9,12,16], the CS decomposition gives a recursive method for 
synthesizing 2-valued and 3-valued n qudit quantum logic gates. In the 2- 
valued case the CSD of a 2" x 2" unitary matrix W reduces to the form 



W 




(3.0.1) 



with each block matrix in the decomposition of size 2" x 2 



in— 1 



In terms of synthesis, the block diagonal matrices in (3.0.1) are quantum mul- 
tiplexers [12]. A quantum multiplexer is a gate acting on n qubits of which 
one is designated as the control qubit. If the control qubit is the highest order 
qubit, the multiplexer matrix is block diagonal. Depending on whether the 
control qubit carries |0) or |1), the gate then performs either the top left block 
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Fig. 1. 2-valued Quantum Multiplexer M controlling the lower (n — 1) qubits by the top qubit. The slash 
symbol (/) represents (n — 1) qubits on the second wire. The gates labeled +1 are shifters (inverters in 
2-valued logic), increasing the value of the qubit by 1 mod 2 thereby allowing for control by the highest 
qubit value. Depending on the value of the top qubit, one of Ut is applied to the lower qubits for t 6 {0, 1}. 

or the bottom right block of the nxn block diagonal matrix on the remaining 
{n — 1) qubits. A circuit diagram for a n qubit quantum multiplexer with the 
highest order control qubit is given in figure 1. Observe that we decomposed 
arbitrary quantum multiplexer to single qubit gates and n qubit standard con- 
trolled gates. The controlled gates execute the operator in the box when the 
controlling qubit has values 1 (mod 2). Such a quantum multiplexer can be 
expressed as 

/ ^0 \ 

(|ai)® las) «)•••«) |an)) (3.0.2) 

where |aj) is the i-th qubit in the circuit, and both block matrices Uq and Ui 
are of size 2""^ x 2"~^. Depending on whether \ai) — |0) or |ai) — |1), the 
expression (3.0.2) reduces to 



|0) ® Uq (laa) las) «)•••«) \an)) 



(3.0.3) 



or 



respectively. 



Ui{\a2) 



03) 



® \an)) 



(3.0.4) 



The cosine-sine matrix in (3.0.1) is realized as a uniformly (n — 1)- controlled 
Ry rotation gate, a variation of the quantum multiplexer. As shown in figure 2, 
a uniformly (n — 1) -controlled Ry rotation gate Ry is composed of a sequence of 
(n — l)-fold controlled gates all acting on the highest order qubit, where 



Rl 



cos 9i — sin 9i 
sin 9j cos 9i 



(3.0.5) 



The control selecting the angle 6i in the gate Ry^ depends on which of the 
(n — 1) basis state configurations the control qubits are in at that particular 
stage in the circuit. In figure 2, the open controls represent the basis state |0) 
and a filled in control represents basis state The i-th {n — l)-controlled 
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Fig. 2. A uniformly (n — l)-controllcd _Ry rotation for 2- valued quantum logic. The lower (n — 1) qubits 
are the control qubits represented on the left hand side by the symbol / on the second wire. The o control 
turns on for control value |0) and the • control turns on for control value |1). It requires 2'^~^ one qubit 
controlled gates iiy* to implement a uniformly (n — l)-controlled Ry rotation. 
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Fig. 3. A control by input value (mod 2) realized in terms of control by the highest value 1 (mod 2). 



Fig. 4. A uniformly 2-controlled Ry rotation in 2-valued logic: the lower two qubits are the control qubits, 
and the top bit is the target bit. 

gate Ry^ may be expressed as 
/ 



cos^i — sin^j 
sin Oj cos 9i 



® (|a2) ® ■■ ■® |a„)) 



(3.0.6) 



with 9i taking on values from the set {^qj ^Ij • • • i ^2"-i-i} depending on the 
configuration of (|a2) ® ■ ■ ■ ® \cin))-i resulting in a specific Ry^ for each i. 

As an example, consider the 3 qubit uniformly 2-controlled Ry gate controlling 
the top qubit from figure 4. Then the action of Ry^ on the circuit is 



cos^j — sin^i 
sin 9i cos 9i 



\ai) ® (|a2) ® las)) 



(3.0.7) 



with 9i G {^0; 9i, 92, 9^}. As \a2) ® {a^) takes on the values from the set 
{|0) ® |0) , |0) (g) |1) , |1) ® |0) , |1) ® |1)} in order, the expression in (3.0.7) re- 
duces to the following 4 expressions respectively. 



cos^o — sin^o 
sin 9o cos 9o 



|ai)®(|0)® |0)) 



(3.0.8) 
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Fig. 5. 3-valued Quantum Multiplexer M controlling the lower (n— 1) qutrits via the top qutrit. The slash 
symbol (/) represents (n — 1) qutrits on the second wire. The gates labeled +2 are shift gates, increasing 
the value of the qutrit by 2 mod 3, and the control turns on for input |2}. Depending on the value of the 
top qutrit, one of Ut is applied to the lower qutrits for t € {0, 1, 2}. 



'^COS^l — sin^i 

I 

sin 01 cos 9i 



V 

^ cos 62 — sin^^2 

I 

sin 62 cos 62 



®m®\i)) 
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(3.0.9) 



(3.0.10) 



V 



cos 9-i — sin 9-1, 
sin 6*3 cos 6*3 



|ai)®(|l)® |1)) 



(3.0.11) 



4 CSD Synthesis of 3-valued (ternary) Quantum Logic Circuits 



In the 3-valued case, two applications of the CSD are needed to decompose a 
3" X 3" unitary matrix W to the point where every block in the decomposition 
has size ?>^~^ x 3""^ [8]. Choose the parameters m and r given in (2.0.1) as 
m = 3" and r = 3"-\ so that m - r = 3" - 3"-^ = 3"-^(3 - 1) = 3""^ • 2. 
The CS decomposition of W will now take the form in (2.0.2), with the matrix 
blocks U and X of size 3""^ x 3""^ and blocks V and Y of size 3"~^ - 2 x 3""^ -2. 
Repeating the partitioning process for the blocks V and Y with m = 3""'^ ■ 2 
and r = 3"~^, and decomposing them with CSD followed by some matrix 
factoring will give rise to a decomposition of W involving unitary blocks each 
of size 3"~^ as follows. 
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(4.0.1) 
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Fig. 6. A uniformly (n — l)-controlled Rx rotation. The lower (n — 1) qutrits are the control qutrits 

represented on the left hand side by the symbol / on the second wire. The controls o, and o turn on for 
inputs |0), |1), and |2) respectively. It requires 3"~^ one qutrit controlled gates to implement a uniformly 
(n — l)-controlled Rx or Rz rotation. 



f2 — ^ 2 — 



Fig. 7. A control by the value (mod ,3) realized in terms of control by the highest value 2 (mod 3). 



M — 1 



Fig. 8. A control by the value 1 (mod 3) realized in terms of control by the highest value 2 (mod 3). 

with 
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(4.0.2) 
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We realize each block diagonal matrix in (4.0.2) and (4.0.3) as a 3-valued 
quantum multiplexer acting on n qutrits of which the highest order qutrit 
is designated as the control qutrit. Depending on which of the values |0), 
|1), or 1 2) the control qutrit carries, the gate then performs either the top 
left block, the middle block, or the bottom right block respectively on the 
remaining n — 1 qutrits. Figure 5 gives the layout for a n qutrit quantum 
multiplexer realized in terms of Muthukrishnan- Stroud (MS) gates. The MS 
gate is a d-valued generalization of the controUed-not (CNOT) gate from 2- 
valued quantum logic, and allows for control of one qudit by the other via 
the highest value of a d- valued quantum system, which in the 3-valued case is 
2 [10]. 
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Fig. 9. An n qudit hybrid quantum multiplexer, here realized in terms of Muthukrishnan-Stroud (d- valued 
controlled) gates. The top qudit has dimension dt and controls the remaining (n — 1) qudits of possibly 
distinct dimensions which are represented here by the symbol (/). The control turns on for input value 
\di ~ 1} mod di of the controlling signal coming from the top qudit. The gates +{di — 1) shift the values of 
control qudit by {di — 1) mod d,. 



The cosine-sine matrices are realized as the uniformly (n — l)-controlled 
and Rz rotations in R^. Similar to the 2- valued case, each R^ and R^ rotation 
is composed of a sequence of {n — l)-fold controlled gates R^' or Rp, where 



Rl^ 



1 

COS 9i — sin 9, 
sin^,- cos^, 



^ cos 0i — sin 0i ^ 



Rt 



sm (, 




cos 0j 
1 



(4.0.4) 



Each R^' or Rf^ operator is applied to the top most qutrit, with the value of 
the angles 9i and (pi determined by the (n — 1) basis state configurations of the 
control qutrits. A uniformly controlled Rx gate is shown in figure 6. Figures 7 
and 8 explain the method to create controls of maximum value. The value of 
the control qubit is always restored in figures 7 and 8. 



5 Synthesis of Hybrid and d-valued Quantum Logic Circuits 



It is evident from the 2 and 3-valued cases above that the CSD method of 
synthesis is of a general nature and can be extended to synthesis of d-valued 
gates acting on n qudits. In fact, it can be generalized for synthesis of hybrid 
n qudit gates. We propose that a (^1^2 ■ ■ ■ dn) x (^1^2 ■ ■ ■ dn) block diagonal 
unitary matrix be regarded as a quantum multiplexer for an n qiidit hybrid 
quantum state space 7i — Hd^ <8) Hd^ <S> - ■ ■ <^ H-dn ? where Ti^. is the state space 
of the i qudit. 

Moreover, consider a cosine-sine matrix of size ((^1^2 ■ ■ ■ dn) x {did2 ■ ■ ■ dn) of 
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Fig. 10. A hybrid uniformly (n — l)-controlled Givens rotation. The lower (n — 1) qudits of dimensions 

d2, d:j, . . . , di-i, rfi+i , . . . ,dn, respectively, are the control qudits, and the top is the target qudit of dimension 

( k) 

di. The control gate turns on whenever the control qudit of dimension d; takes on the value k (mod) 
di. 



the form 



hp ^ 
C -S 
S C 



7, 



(5.0.1) 



1/ 



with Ip and Ig both some appropriate sized identity matrices, C — diag(cos^i, 
cos 02, ■ ■ ■ , cos 9t) and 5" = diag(sin 9i, sin $2, . . . , sin such that sin^ 9^ + 
cos^ 6i = 1 for some with 1 < i < t, and p + g + 2t = (^1^2 • • • c?„). We 
regard this matrix as a uniformly controlled Givens rotation matrix, a gener- 
ahzation of the Ry, Rx, and Rz rotations of the 2 and 3- valued cases. A Givens 
rotation matrix has the general form 



1 ... 



... 



. . . cos ^ 



sin ^ . . . 



. . . sin ^ . . . cos 9 



... 







(5.0.2) 



7 



where the cosine and sine values reside in the intersection of the i-th and 
j-th rows and columns, and all other diagonal entries are 1 [7]. Hence, a 
Givens rotation matrix corresponds to a rotation by some angle 9 in the ij'-th 
hyperplane. 
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Fig. 11. A control by the value k (mod rfj) realized via control by the highest value (d; 



— 1) (mod di). 



Based on the preceding discussion, we give in theorem 5.1.1 below an iterative 
CSD method for synthesizing a n qudit hybrid quantum circuit by decom- 
posing the corresponding unitary matrix of size {did2 . . . dn) x {did2 ■ ■ ■ dn) in 
terms of quantum muhiplexers and uniformly controlled Givcns rotations. As 
a consequence of this theorem, we give in lemma 5.1.1 a CSD synthesis of 
a quantum quantum logic circuit with corresponding unitary matrix of size 
X d". The synthesis methods given above for 2- valued and 3- valued circuits 
may then be treated as special cases of the former. 

5.1 Hybrid Quantum Logic Circuits 

Consider a hybrid quantum state space of a n qudits, 7i = 'Hdi®'Hd2®' ' ■^'Hd„, 
where each qudit may be of distinct (i- valued dimension di, < i < n. Since 
a qudit in 7i is a column vector of length did2- ■ - dn, a quantum logic gate 
acting on such a vector is a {did^ ■ ■ ■ d^v) x {did2 ■ ■ ■ dn) unitary matrix W. 
We will decompose W, using CSD iteratively, from the level of n qudits to 
{n — 1) qudits in terms of quantum multiplexers and uniformly controlled 
Givens rotations. However, since the d-valued dimension may be different for 
each qudit, the block matrices resulting from the CS decomposition may not 
be of the form x for some d. Therefore, we proceed by choosing one of 
the qudits, q. of dimension di, to be the control qudit and order of the basis of 
Ti in such a way that is the highest order qudit. We will decompose W with 
respect to q. so that the resulting quantum multiplexers are controlled by 
and the uniformly controlled Givens rotations control q. via the remaining 
{n — 1) qudits. We give the synthesis method in theorem 5.1.1. 

Theorem 5.1.1: Let W he an M x M unitary matrix, with M = did2 . . . dn, 
acting as a quantum logic gate on a quantum hybrid state space Ti = Tidi ® 
'Hd2 <8) • • • <8) Ti.d„ of n qudits. Then W can he synthesized with respect to a 
control qudit c^^ of dimension di, having the highest order in Ti, iteratively 
from level n to level {n — 1) in terms of quantum multiplexers and uniformly 
controlled Givens rotations. 

Proof. 

Step 1. At level n, identify a control qudit q. of dimension di. Reorder the 
basis of Ti so that q. is the highest order qudit and the new state space 
isomorphic to is = Ti^. ® Hd2 <8) • • • (8) Hd^ <8) • • • (8) Hd^- 
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If we choose values for the CSD parameters m and r as m = {did2 ■ ■ ■ d,!) 
and r = {d,id2 . . . di^idi^i . . . dn), then m — r = di . . . dj^idj+i . . . dn{di — 1). 
Decomposing W by CSD, we get the form in (2.0.2) with the matrix blocks U 
and X of size r x r and blocks V and Y of size {m — r)x{m — r). Should m — r 
not have the factor (di — 1), we would achieve the desired decomposition of W 
from level of n qudits to the level of (n — 1) qudits in terms of block matrices 
of size r X r. The task therefore is to divide out the factor (rfj — 1) from m — r 
by an iterative lateral decomposition described below, that uses the CSD to 
cancel {di — 1) from m — r at each iteration level leaving only blocks of size 
r X r. 

For step 2 of the proof below, we will say that a matrix with k rows and k 
columns has size k instead of k x k. 

Step 2. Iterative Lateral Decomposition: For the unitary matrix W of size 
M, we define the j-th lateral decomposition of as the CS decomposition of 
all block matrices of size other than r that result from the (j — l)-st lateral 
decomposition of W: 

ForO <j< {di-2), set 
mo = {did2 ...dn) 
tq = {did2 . . . di-idi+i ...dn) 
If J = 

Apply CSD to W 
Else set 

ruj = mo - j • ro 
rj = ro 

rrij - rj = mo - (j + l)ro 

={did2 . . . di^idi+i . . . dn) [di - {j + 1)] 

mj - 2rj = mo - (j + 2)ro 

=(^1^2 ■ ■ ■ di-idi+i ...dn) [di - {j + 2)] 

Apply CSD to matrix blocks of size other than ro from step j — 1 
End If 
End For. 

When j — 0, we call the resulting 0-th lateral decomposition the global decom- 
position. Note that if di — 2, then the algorithm for the lateral decomposition 
stops after the global decomposition. This suggests that whenever feasible, the 
control system in the quantum circuit should be 2- valued so as to reduce the 
number of iterations . Below we give a matrix description of the algorithm. 

For j — 0, the 0-th lateral decomposition of W will just be the CS decompo- 
sition of W. 

(5.1.1) 
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where 



4(0) _ 
^0 — 
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q(0) ^(0) 
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(0) 











with U^^\ Xf, Cf\ 



and S'o°'' aU of the desired size ro, while 



(0) 



and 



(0) 



are of size mo — tq. The superscripts label the iteration step, in this case 
j' = 0. The subscript is used to distinguish between the various matrix blocks 
[/, V, X, Y, C, S, that occur at the various levels of iteration. The 0-th lateral 
decomposition in the form from equation (5.1.1) is called the global decompo- 
sition of W. 



For 



1, we perform lateral decomposition on the blocks Vq^^ and Yq^^ of the 



block matrices Aq' and Cq"' respectively, the only blocks of size other than tq 



and 

resulting from the 0-th lateral decomposition given in (5.1.1). In both cases, 
set nil = ^0 — ^0 and ri = ro so that mi — ri = mo — 2ro. For Vq^^ this gives 
the decomposition 



/MO) 
^0 
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/ ^(1) q(1) 
'-'0 ~'-'o 
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(1) 







(5.1.2) 



with Uq'\Xq'\Cq'^ and Sq^' all of size ro, and Vq^' and Fo''^'' of size mi — ri. 
All three matrices residing in the lower block diagonal of the matrix (5.1.2) 
are the same size. Therefore, by introducing identity matrices of size ro and 
factoring out at the matrix block level, Aq'^ will be updated to 



(1) 



A 



(0) 



-^0 '-"0 



(5.1.3) 



where 
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A similar lateral decomposition of the block Y^^ will update C^^ in (5.1.1) 
to 



Cf = A\"B\"C[ 



(1) nWr^W 



(5.1.4) 
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where 
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cr -s\ 
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r/^^ y 



For iteration j ^ 0, perform lateral decomposition on the total 2^ blocks 
Vk~^\ Yk~^\ where < A; < (j - 1), that occur in the global decomposition 
at the end of iteration (j — 1). For each V^^~^\ Y^^~^\ set Vj — Tq, rrij — 
ruj^i — Tj-i — mo — {j + l)ro. For each V^^~^\ the lateral decomposition at 
level j will give the following 
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(5.1.5) 

where the A^-' is the block diagonal matrix of size of j ■ tq arising from the 
lateral decomposition in the previous (j - 1) steps. The blocks f/f.f\ , d^) 
and Sf) arc all of size tq, for < A;' < j. The blocks V^^'' and are of size 
mj — Tj. The three matrices residing in the lower block diagonal of the matrix 
(5.1.5) are all of same size. Therefore, by introducing identity matrices of size 
j • To and factoring out at the block level, will be updated to 



^fe' — ^fe' ^k' ^k' 



where 



^k' 



A(^-i) 
U^,^ 




^k' 



R(i) 











\ 





^k' 


qU) 








qij) 







V 








Imo-{j+2)ro J 



^k' 



Ij-ro 

x;^,^ 

Y^P 



For the next iteration, set k — k' and iterate. Upon completion of the lateral 

decomposition, repeat steps 1 and 2 for the synthesis of the circuit for the 
remaining (n — 1) qudits, with the restriction that each gate in the remaining 
circuit be decomposed with respect to the same control qudit identified in step 
1. 
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Since the basis for Ti was reordered in the beginning so that the control qudit 
was of the highest order, the block diagonal matrices with all blocks of size 
Tq X To are interpreted as quantum multiplexers and the cosine-sine matrices 
are interpreted as uniformly controlled Givens rotations. In figures 9 and 10, 
we present the circuit diagrams of a hybrid quantum multiplexer and a uni- 
formly controlled Givens rotation, respectively. A uniformly controlled Givens 
rotation matrix on n qudits can be realized as the composition of various 
{n — l)-fold controlled Givens rotation matrices, G^f^j, acting on the top most 
qudit of the circuit with the angle of rotation depending on the basis state 
configuration, in their respective dimensions, of the lower {n — 1) qudits. 

5.2 d-valued Quantum Logic Circuits 

Given the hybrid n qudit synthesis, the case of rf-valued synthesis becomes 
a special case of the former since by setting all di = d, the state space 7i ~ 
'Hd]_®'Hd2®- ■ -^Ti-du reduces to the state space TY®". Unitary operators acting 
on the states in V,®'^ are unitary matrices of size d^xd^. We give the following 
result for o?-valued synthesis. 

Lemma 5.1.1: A d-valued n qudit quantum logic gate can he synthesized in 
terms of quantum multiplexers and uniformly controlled Givens rotations. 

Proof: Since all the qudits are of the same dimension, there is no need to 
choose a control qudit. In the proof of theorem 5.1.1, set d^ = d for all i. Then 
M — did2 . . .dn — d^. For iteration j = of the lateral decomposition, set 
mo = d", ro = d""-^, so that mo - ro = d'^'^id - 1). For < j < (d - 2), set 
rj — rQ — d^~^, and rrij — rrij-i — rj-i — d^'^ld — {j + 1)). 

For the d-valued case, we note that there are a total of d"'~^{2^'^ — 1) one 
qudit Givens rotations in the circuit at the (n — 1) level, each arising from 
the Z^iio^'' 2* = 2'^~^ — 1 uniformly controlled Givens rotations in the CS 
decomposition of an n qudit gate. Moreover, in each uniformly controlled 
Givens rotation, there are (n — l)d"'~^ control symbols of which (n — l)(i"~^ 
correspond to control by the highest value of d — 1. The latter controls do not 
require shift gates around them to increase the value of the signal qudit to d—1. 
Hence, there arc (n— — (n— l)(i'*~^ = {n—l){d^~^—d"'~'^) control symbols 
that correspond to control by values other than d—1 and therefore need two 
shift gates (fig. 11) around them. This gives the total number of one qudit 
shift gates in each uniformly controlled rotation to be 2(n — l){d'^~^ — o?'*"^), 
whereby the total number of one qudit shifts and Givens rotations in the 
circuit at the (n-1) level is 2(n - l)(V-i -d"-2)(2<^-^ - 1) + d"-i(2''-i - 1) = 
(2'i-i - 1) [2{n - l)(ci"-i - d""-^) + 

There are 2'^"^ quantum multiplexers in the decomposition, each consisting of 
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a total of 2d shift and controlled gates. Hence, there are a total oi d ■ 2'^ one 
qudit and controlled gates in the (n — 1) level circuit. This gives a total, worst 
case, one qudit and controlled gate count in the circuit at level (n — 1) to be 
- 1) [2{n - - + (^"-1] + d-2'^. 



6 Conclusion 

We have shown that the method of CS decomposition of unitary matrices used 
in 2-valued and 3-valued quantum logic synthesis is a special case of a general 
synthesis method based on the CSD. We give an algorithm for this general 
method that allows us to synthesize n qudit hybrid and d- valued quantum logic 
circuits in terms of quantum multiplexers and uniformly controlled Givens 
rotations. 
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